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Abstract 

Breaking of an atomic chain under stress is a collective many-particle tun- 
neling phenomenon. We study classical dynamics in imaginary time by using 
conformal mapping technique, and derive an analytic formula for the prob- 
ability of breaking. The result covers a broad temperature interval and in- 
terpolates between two regimes: tunneling and thermal activation. Also, we 
consider the breaking induced by an ultrasonic wave propagating in the chain, 
and propose to observe it in an STM experiment. 

PACS numbers: 05.30.-d, 61.16.Di, 11.17.+y 



1 



Typeset using REVTgX 



Introduction. Among various examples of quantum decay in condensed matter |lj of special 
interest are those where the tunneling system is strongly coupled to a macroscopic system. 
Several such systems were studied, including metastable current states of Josephson junc- 
tions P|, transport in tunnel junctions [||], quantum diffusion and properties of two-level 
systems coupled to thermal bath ||]. Our purpose here is to discuss yet another system 
where decay is a many-particle effect: breaking of elastic string stretched by an external 
force. Our central example will be a polymer chain teared under stress, however the for- 
malism we present is more general. For instance, it describes critical velocity of a superffuid 
flow in a long channel and breakdown of a supercurrent carried by a thin wire or by many 
small capacitance Josephson junctions connected in series. 

For the string, the leading contribution to the tunneling exponent is due to a collective 
motion during breaking |5|]. According to the Caldeira-Leggett theory ||, since the string 
corresponds to their "ohmic" case and is described by a harmonic Hamiltonian, the internal 
degrees of freedom can be integrated out to produce an effective action for which a classical 
trajectory in imaginary time has to be determined. Although in principle there is no difficulty 
in carrying out this program, we chose a different method. The "collective motion" is just 
a propagation of acoustic waves, and it can be treated exactly: in imaginary time the 
wave equation becomes two dimensional Laplace's equation that can be solved by conformal 
transformations. The conformal mapping method gives a clear intuitive picture of the entire 
system motion in imaginary time and allows to treat in parallel several related problems: 
tunneling at zero temperature and at finite temperature, and also tunneling caused by 
an alternating tension. The latter problem is of a particular practical interest because it 
corresponds to the breaking induced by an ultrasonic wave propagating in a polymer. 

We consider the limit of small tension, when the tunneling results from coherent mo- 
tion of many atoms. A simple argument that allows to estimate the number of moving 
atoms and the exponent in the breaking probability was given by Dyakonov ||. During 
the time of tunneling t, acoustic waves propagate at the distance ct, where c is the sound 
velocity. At small tension, the time t is big, and thus the mass involved in the tunneling, 
m* = pet, is much bigger than the single atom mass. (Here p is the density of the chain.) 
One can then describe dynamics approximately as tunneling of the mass m* under a tri- 
angular barrier U(x) = E — fx, where E is the energy of the broken bond and / is the 
tension. Then it is straightforward to estimate the action S ~ f~~ l rrhj 2 E^ 2 and the time 
t ~ f^rruj* ''E^ 2 '. By solving the latter equation together with the selfconsistency relation 
m* = pet, one finds t ~ pcE /f 2 , ~ p 2 c 2 E / f 2 , which gives the probability of breaking 
W ~ exp(— constE^ pc/hf 2 ). The argument of Dyakonov shows that collective effects are 
essential, as grows at small tension. In this form, however, it only gives a rough estimate 
of lnVF, leaving an unknown constant in the exponent to be determined by a more precise 
calculation. 

Discussion of experiment. Beyond a theoretical interest, the quantum problem of string 
breaking is relevant for the physics of polymers, since at low temperature this process de- 
termines maximal stress o max that a polymer material can sustain. At high temperature 
the breaking is thermally assisted: o~ max (T) ~ expEo/T. At low temperature, a saturation 
of the raise of a max iT) was reported however, the relation with tunneling has not been 
clarified. It would be of interest to check whether the saturation indicates tunneling by 
using, e.g., the conventional technique of isotopic substitution. Since \iaW scales as l/pc, 
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there is an isotope effect in the quantum regime: In W ~ VM, which disappears in the 
thermally assisted regime. Also, one could look at breaking enhancement in presence of an 
alternating stress caused, e.g., by a microwave radiation. In this case, the characteristic 
feature to observe is the threshold frequency Qq where the breaking enhancement described 
by our Eq. (|T9"D begins. The threshold is set by inverse tunneling time t' 1 ~ f 2 /pcEo, and 
thus it is a function of the tension / and can be much lower than Debye's frequency uu, the 
characteristic threshold for the frequency effect in the thermally assisted regime. 

Also, in the light of recent progress in manipulating single atoms by a scanning tunneling 
microscope 0, it is of interest to discuss an STM experiment with a single polymer molecule 
suspended between the surface and the STM tip acting as a probe. In this setup one 
can study the molecule breaking at small and at large temperature, where our calculation 
predicts different regimes. In addition, there is a possibility to combine constant force with 
an AC tension field caused, e.g., by an AC signal on the tip, and to look at resonance effects 
related with exciting a standing ultrasonic wave in the molecule. Due to a large acoustic 
mismatch with the surface, reflection of the ultrasonic wave at the ends of the molecule is 
almost perfect. Thus the molecule of length L works like a high quality phonon resonator 
with resonance frequencies ui n = nnc/L corresponding to sinusoidal standing waves. For 
example, in a chain of 1000 atoms one estimates u)\ as 10~ 3 co>£> « 5 — 10 GHz. A signal in 
this frequency range presumably can be supplied by coupling the STM tip to a microwave 
source. 

Breaking then will occur at the nodes of the wave, where local tension is maximal, and 
this will result in simple ratios of broken halves: 1:1, 2:1, 3:2, etc. On the other hand, at high 
frequency there is no particular spatial structure of the breaking. Like in the bulk polymer, 
the experiment is facilitated by anomalously low frequency threshold Qq ~ t^ 1 above which 
the breaking is exponentially enhanced. 

Result and method. The string is described by the wave equation that in imaginary time 
turns into Laplace's equation. We solve the latter by a conformal mapping and derive a 
formula for the probability of breaking 

*-A-p(-3f). 

where A ~ UDe~ E °^ D is the probability of a single atom tunneling at one atomic spacing. 
Let us note an agreement with the qualitative estimate given before. We generalize this 
result to a finite temperature and derive an exact formula that covers a broad temperature 
range, including thermal activation at high temperatures. 

We shall use the following picture of tunneling. The energy U (x) of a bond falls off as 
its length x exceeds several angstroms. Therefore, due to zero-point fluctuations, the bond 
occasionally makes excursions in the state where it is virtually broken, with the frequency 
~ A. If there were no tension, after every such excursion the bond would return back to its 
equilibrium length. However, in the presence of the force stretching the chain, after the bond 
is virtually broken, the ends of the tear can start moving apart in the classically forbidden 
region, until the distance between them reaches ~ Eq/ f, where the barrier ends and the bond 
is really broken. This second stage takes much longer than A -1 as it involves motion of many 
atoms at large distance. Thus, we have two stages of tunneling: (1) fast virtual disappearance 
of one bond; (2) slow collective motion under applied force traversing classically forbidden 
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region. After we use this picture to do the calculation, we check selfconsistency of the 
assumptions, and confirm validity of the approach for the tension small on the atomic scale. 

Clearly, since the first stage is essentially one-particle and does not depend on the tension, 
it only contributes a prefactor ~ A 1 / 2 of the exponentially small amplitude due to the second 
stage. The latter can therefore be analyzed using the method of sudden switching ||. Let 
us consider a Hamiltonian 



k=l...N 



+ U{ qi ) + U{L-q N ) 



where M is atomic mass, U (x) is the interaction of neighbors, are atoms' coordinates. The 
ends of the chain are fixed at x = 0, L, which determines interatomic spacing a = Lj (N + 1) 
and the tension / = U'(a). Our / is small, so a is close to the equilibrium spacing a for 
which U'(ao) = 0, and we can use harmonic approximation 

U (x) = -Eo + ^muj 2 {x - a f (2) 

for all bonds except the one that breaks. For the breaking bond, we practically do not need 
exact profile of the potential U(x), as it is important only during a very short time interval 
~ A -1 at the beginning of the tunneling. Therefore, as we discussed, it contributes only to 
the prefactor of Eq.(|). According to |§, we write the rate of transition to the broken state 
as 

oo 

W = A J (exp(~H't))dt, (3) 

— oo 

where the Hamiltonian ri' describes the chain with one bond removed: TC' — TC — U (xj — 
Xj-i) — E . (We subtract E to assure energy conservation.) The average in Eq.(||) is taken 
over the initial state, i.e., the ground state of H. 



Following the usual treatment ||10|| , we consider the tunneling as a motion in imaginary 
time: first, we evaluate the average in (^) at imaginary values of t, then continue it from the 
upper half-plane down to the real axis and do the integral over t by the saddle point method. 
To compute the average, we write exp(—fH'/fi) averaged over the initial wave function as 
the path integral 

(exp(-~#)) = e- E ^ n J V[q] exp (~S it/2 [q\) (4) 

Here 5j r is the classical action of the chain with one bond broken during the time interval 
\—zt, it]: 



S iT [q] = J [n- 6(t 2 - t 2 )U( gj - q^)\ dt , (5) 

where j is the number of the broken bond. The integral (Q) is taken over paths going from 
t = — oo to t — oo, i.e., over all states of the chain. Since at small tension the 

tunneling involves only long wavelength distortions, we go to the continual limit: 
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P-2 . P c2 2 , 

2 U + ^ u *~f Ul 



dxdt , (6) 



where x{k) = a>ok, and the atoms' coordinates qk are written through the displacements 
relative to the equilibrium positions: qk(t) = u(x,t) + a k. Together with the boundary 
condition corresponding to the broken bond at x — 0, this defines a linear problem. We 
solve it and evaluate the functional integral by the saddle point method. 

Beginning from here we use the units % = p = c = 1, and recover the usual units only 
when necessary. Classical equation of motion obtained from Eq.(^|) is 

d 2 u d 2 u 

Ih 2 + ~W ' { } 

We assume that the breaking occurs far from the ends and treat the chain as infinite. Then 
the Laplace's equation is supplied with the boundary conditions: 

(i) u(x,t) — > u (x,t) = fx at \x\, \t\ ^> r, i.e., the chain is uniformly stretched away from 
the tear; 

(ii) du/dx = at x = ±0, — r < t < r, i.e., the ends of the tear at x = are free during the 
tunneling. 

By analogy with 2D electrostatics, the boundary value problem can be solved by a proper 
conformal mapping. For that, let us introduce complex variable z = x + it and consider 
analytic function 



w 



(z) = (z + Vz 2 + t^)/t . (8) 



The cut [— 2T, 2t] in the z— plane is mapped on the unit circle \w\ = 1, and thus we come to 
Laplace's equation in the domain \w\ > 1 with the boundary conditions 

du 

(i) u = uq(w) at \w\ 3> 1; (ii) —, — , = at \w\ = 1 (9) 

o\w\ 



u (z) = i/rRe (w + — ) = fReVz 2 + r 2 . (10) 



It has a solution 



2" V w 

Due to conformal invariance, the action can be evaluated in the w— plane: Si T [u — uo] = 
—iif 2 T 2 /2. Then we continue Si T to the real axis, substitute it in Eq.(^), do the integral 

W ~ J drexp [^-\f^ 2 ~ 2iE rj (11) 

and get Eq.(|l]) for the probability. 

With the result flT0D, we can check the selfconsistency of the method. During the tun- 
neling, the width of the tear u(x = +0,t) — u(x = — 0,t) = ^:Vt 2 — t 2 is estimated as E Q /f, 
as the saddle point evaluation of the integral (0) sets r = 2ipcE /nf 2 . If the tension / is 
small in atomic units, the tear is much wider than the interatomic spacing, except for t very 
close to ±r, and thus indeed it can be described by the free end boundary condition at the 
cut. 
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Finite temperature. Temperature dependence of the breaking begins very early, at the 
temperature set by the inverse tunneling time, T c ~ H/r = hf 2 /pcE . If the tension is weak, 
T c can be much smaller than the usual temperature Tiujd of transition to the thermally 
assisted tunneling. For such temperature the breaking can be studied by a generalization of 
the zero temperature calculation. 

The functional integral treatment of tunneling at finite temperature amounts to inte- 
grating over periodic trajectories with the imaginary period i(3 = i/T [|10|| . For the chain, 
we have to evaluate the path integral Eq.(|]) taken over all functions u(z) periodic in the 
z— plane with the period i/3: u{z + i(3) = u(z). Repeating the steps that lead to Laplace's 
equation, we get a boundary value problem for the function u(z) whose normal derivative 
vanishes near a periodic system of cuts, z E [—it + i/3n, ir + i/3n\. (Condition at z = oo 
remains unchanged.) 

Corresponding conformal mapping can be constructed in two successive steps. First, let 
us consider the function r](z) = exp(2irTz) that transforms the strip < Imz < (3 to the 
whole complex plane, and maps z = oo to two points: r\ = 0, rj = oo. The cut is mapped to 
the arc of the unit circle that goes from a = exp(27riTr) to a counterclockwise. The second 
function we take is 

rj-a 

w(V) = 7 ■ 

arj — 1 

In the w— plane the two images of z = oo are ui\ = a, w 2 = a, and the cut goes along the 
negative real half axis w < 0. Then the boundary value problem is solved by 




with the standard branch of the square root. Computation of the action Si T [u — u Q ] is 
facilitated by comparing with analogous problem of 2D magnetostatics: two opposite charge 
magnetic monopoles situated at the points w 1>2 near a superconducting slit going along the 
real axis, from to the left. In this formulation, the action corresponds to the interaction 
energy of the monopoles and their images. Straightforward calculation gives 

S iT [u — Uq] = Alncos 2 (7rTr) , A 



2nT 2 



To obtain the probability W, we take ir back to the real axis, substitute S T in Eq.(|), and 
do the integral: 

W ~ f e2lE0TdT 22A r ^ + ) r (A-) 

J cosh 2A (7rTr) 2vrT T(2A) ' K } 

A± = A ± iEo/nT. This expression covers a large temperature range, up to T ~ up. By 
using Stirling formula we rewrite the result: 

lniy(/,T) = -^^2tan-^-iln(e 2 + l)) , 

where £ = 2TEopc/hf 2 . Asymptotically, at high temperature, £ 3> 1, we get Arrhenius law 
W ~ \e~ E °/ T , and in the opposite limit, £ <C 1, Eq. ([T2|) matches Eq.(P. 
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Breaking due to alternating tension. Now, let us discuss how the breaking can be stimulated 
by a sound wave. Local tension in the wave is varying along the string, and so is the breaking 
probability. To study the spatial dependence, let us put a phase in the wave, 

fo 

Uq(x, t) — — — cos(f2(a:/c — t) + (p) , (13) 

where fo is the tension amplitude, and then study breaking at x = at different values of 
0. In the z— plane 

u (z) = ^Ree^ z -^ , (14) 

where k = Q/c. Let us assume zero temperature, then we have Laplace's equation for 
the displacement field u(z) with the boundary conditions identical to that we had in the 
T = constant force problem. By the conformal transformation Eq. (^) we go to the domain 
\w\ > 1 in the w— plane, z = t(w — l/w)/2, where we get the boundary value problem @ 
with uo(w) corresponding to Eq. (fH]) . Solution is readily obtained by writing Laurent series 
for u(w) and uq(w), solving separately for each term, and then collecting the series: 

2fo 00 2 m +i 

11(10) = ito(tu) + — -^sin0Re ^ —^-J m {iVtT) , (15) 

^ m=l W 

where J m (x) are Bessel functions. The action is 

S iT [u-u ] = -27T-^ sin 2 ^2 m\J m (iQT)\ 2 (16) 

m=l 

This sum can be done in the following two cases. 

I. (low frequency: Qr <^ 1) Here J n (x) = 0(x n ), and only the m = 1 term survives: 

^ = -^sin 2 (17) 
2pc 

This expression is just what one has for constant tension with the local value / = /osin0. 
So, in this limit the breaking occurs predominantly at the nodes of the displacement, where 
the tension is maximal. 

II. (high frequency: f2r 3> 1) Here we use the x 3> m asymptotics, J m (ix) ~ i m e x / \f2nx, 
and get 

pc\l 

So, the probability of breaking in this limit is 



W ~ exp 



E ( E pdl 
In 



hQ \ /q sin^ 



(19) 



Let us note a similarity between Eq.(|T9"D and the ionization rate of an atom in a low frequency 
electric field / cos fit: W ~ f 2n , where n = Ry/hQ is interpreted as the number of absorbed 
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quanta and f 2 as the single quantum absorption rate. This interpretation stands meaningful 
for the string as well. 

To summarize, the frequency dependence of the breaking probability is as follows: at 
Q <C Q = r" 1 = f^/pcEo it is given by (|l|), and the breaking occurs at the nodes of the 
wave. At Q > Qq the frequency dependence begins and the probability sharply increases 
following Eq.(^). The spatial dependence practically vanishes in this limit, since many 
cycles are repeated during the tunneling. 

To conclude, we studied breaking of a polymer chain, and derived an analytic expression 
for the breaking probability that describes transition from quantum tunneling to Arrhenius 
law. Also, we studied breaking induced by an ultrasonic wave, and discussed opportunities 
for experiment with bulk polymers and with an STM used as a probe. 
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